Optimal control of a dengue epidemic model with vaccination* 
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Abstract 

We present a SIR+ASI epidemic model to describe the interaction between human and 
dengue fever mosquito populations. A control strategy in the form of vaccination, to decrease 
the number of infected individuals, is used. An optimal control approach is applied in order to 
find the best way to fight the disease. 
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1 Introduction 

Dengue fever is a vector borne disease, which has become an increasingly public health problem 
that carries a huge financial burden to the governments. Currently, the only way of controlling the 
disease is to minimize the vector population. Dengue vaccine for effective prevention and long term 
control under development, is expected to be the solution. Dengue transcends international borders 
and is emerging rapidly as a consequence of globalization and climate changes. It is a disease of great 
complexity, due to the interactions between humans, mosquitoes, and various virus serotypes as well 
as efficient vector survival strategies. The four serotypes, known as DEN1 to DEN4, constitute a 
complex of flaviviridae transmitted by Aedes mosquitos, specially Aedes Aegypti. Infection by any of 
the four serotypes induces lifelong immunity against reinfection by the same type, but only partial 
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and temporary protection against the others. Sequential infection by different serotypes could lead 
to a more severe dengue episode: dengue hemorrhagic fever (DHF). 

Vector control remains the only available strategy against dengue. Despite integrated vector 
control with community participation, along with active disease surveillance and insecticides, there 
are only a few examples of successful dengue prevention and control on a national scale |T| . To make 
matters worse, the levels of resistance of Aedes Aegypti to insecticides has increased, which imply 
shorter intervals between treatments, and only few insecticide products are available in the market 
because of high costs for development and registration and low returns. 

For long time, the evaluation of global dengue disease burden was limited and the stakeholders 
considered the potential market for the dengue vaccine to be small. By the end of 20th century, with 
the increase in dengue infections as well as the prevalence of all four circulating serotypes, faster 
development of a vaccine became a serious concern [2]- It is agreed that a vaccination program 
not only protects directly the individual, but also indirectly the population, which is called herd 
immunity. As a consequence of vaccination, the occurrence of epidemics would decrease relieving 
health facilities. However, constructing a successful vaccine for dengue has been challenging. Not 
only is the knowledge of disease pathogenesis insufficient, but also the vaccine must protect against 
all serotypes so that the level of DHF doesn't increase. 



2 Optimal control of the epidemiological model 

Two types of population were considered: hosts and vectors. The hosts (humans) are divided into 
three complementary classes: susceptible, Sh(t), individuals who can contract the disease; infected, 
Ih(t), individuals capable of transmitting the disease to others; and resistant, Rh(t), individuals 
who have acquired immunity at time t. The total number of hosts is constant, which means that 
Nh = Sh(t) + Ih(t) + Rh(t)- Similarly, the model has also three compartments for the vectors 
(mosquitoes): A m {t), which represents the aquatic phase of the mosquito (including egg, pupae 
and larvae) and the adult phase of the mosquito, with S m (t) and I m (t), susceptible and infected, 
respectively. It is also assumed that N m = S m (t) + I m (t) . The model is described by an initial value 
problem with a system of six differential equations: 

' = PhNh - (Bf3 mh j^- + fi h + uj S h + cjuRh 
= Bp mh ^-S h - {r, h + n h )I h 

= Vhh + uS h - {<7u + nh)Rh 
' i %r = ^{}- ^ft) (s m + im) - {va + ma) A m 1 1 ! 

^™ = VAAm — (^B f3h m j^- + /x m ) S m 

The recruitment rate of human population is noted by ^hNh- The natural death rate for humans and 
mosquitoes, aquatic and adult phase, is described by the parameters fih, Mm an d Ma, respectively. 
We assume that B is the average daily biting (per day) of the mosquito whereas f3 m h and /3^ m are 
related to the transmission probability (per bite) from infected mosquitoes to humans and vice versa. 
By ip we denote the number of eggs at each deposit per capita (per day) . The recovery rate of the 
human population is denoted by r/h . The maturation rate from larvae to adult (per day) is denoted 
by r\A- The vaccine coverage of the susceptible is represented by u (the control variable). The factor 
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a represents the level of inefhcacy of the vaccine: for a = the vaccine is perfectly effective, while 
(7 = 1 means that the vaccine has no effect at all. The main aim of this work is to study the optimal 
vaccination strategy considering both the costs of treatment of infected individuals and the costs of 
vaccination. So, the objective functional is 



minimize J [u] = / [lilh(t) 2 + 7yu(i) 2 ] dt, 
Jo 



(2) 



where 7/ and jy are positive constants representing the weights of the costs of treatment of infected 
and vaccination, respectively. Let Xi(t), with i = 1, . . . , 6, be the co-state variables. The Hamiltonian 
for the present optimal control problem is given by 



H = Ai 



VhNh - [ Bf3 mh — + fih+u) Sh + cruRh 



+ A 2 



B/3 m h—Sh - {rih + M/0 Ih 



+ A 3 [rjhlh + uS h - (<ru + fi h ) R h ] + A 4 



<P 1 



A, 



kN h 



(S m + Im) - {r)A + MA) A r , 



(3) 



A, 



VAAm [ B(3h m — (- [l m J S n 



+ A 6 



+ irfh + ivu 2 . 



By the Pontryagin maximum principle [3J, the optimal control u* should be the one that minimizes, 
at each instant t, the Hamiltonian given by Q, that is, 

H(x*(t),X*(t),u*(t)) = min H (x*(t), X*(t), u) . 
«e[o,i] 

The optimal control, derived from the stationary condition ^ = and considering < u < 1 , is 
given by 

. L f n (Ai - A 3 ) (S h - aR h ) 
u = mm < 1, max < 0, 



271/ 

Substituting the optimal control u* into the state system (fTJ) and the adjoint system X { (t) = — §§-, 
i.e., 

(Ai - A 2 ) (B(3 mh !ff^j + Ai///, + (Ai - X 3 )u 
-2~iil h + X 2 {rj h + fi h ) - X 3 rj h + (A 5 - A 6 ) (^B/3 hm j^^j 
-Xi<ju + X 3 ([i h + au) 
At^ g feivf m + X 4 (t]a + Ma) - A 5 ?/a 



dt 

d\ 2 

dt 
d\ 3 

dt 
d\j 

dt 
d\s 

dt 

dX 6 
dt 



-X±cp (l - ^) + (A 5 - X 6 )B/3 hm jfc + X 5 fi m 
(Ai - A 2 ) (Bp mh j^j - Xiip (l - ^) + A 6j um, 



we obtain the corresponding x* and A*, i = 1, ... 6, with the help of the transversality conditions 
\*(t f ) = 0, i = 1,...,6 (see [3] for details). 

3 Numerical simulation and discussion 



The simulations were carried out using the following values: Nh = 480000, B = 0.5, /3 m h — 0.3, 
/3 hm = 0.3, (i h = 1/(71 x 365), rj h = 1/3, fi m = 1/10, k = 3, m = 3, N m — m x Nh, ip — 6, and 
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Method 


optimal control 


no control (u = 0) 


upper control (u = 1) 


Direct (DOTcvpSB) 


0.146675 


0.674555 


364.940488 


Indirect (backward-forward) 


0.113137 


0.357285 


365.00046 



Table 1: Values of the cost functional ^ 



tf = 365 days. It was considered that the vaccine is imperfect with a level of inefficacy of a = 0.15. 
The initial conditions for the ordinary differential system were: S/,(0) = 2V/, — 216, i/,(0) = 216, 
Rh(0) = 0, A m — k*Nh, S m (0) — N m and I m (0) — 0. The optimal control problem was solved using 
two methods: direct and indirect. For an introduction to direct and indirect methods in optimal 
control we refer the reader to [4j [5] . The direct method uses the optimal functional fl2J) and the state 
system (JXJ) and was solved by DOTcvpSB 6 . It is a toolbox implemented in MatLab, which uses 
an ensemble of numerical methods for solving continuous and mixed-integer dynamic optimization 
problems. The indirect method we used is an iterative method with a Runge-Kutta scheme, solved 
through ode45 of MatLab. The state system with an initial guess is solved forward in time and 
then the adjoint system with the transversality conditions is solved backward in time. The controls 
are updated at the end of each iteration (see [7] for more details). Figure [T] shows the optimal 
control obtained by the two different approaches. They both seem to have the same behavior. 
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Figure 1: Optimal control with direct and indirect Figure 2: Infected humans using different levels of 
approaches. control. 

Table [1] shows the costs obtained by the two methods, in three situations: optimal control, no 
control (u(t) = 0) and upper control (u(t) = 1). Figure [5] shows the number of infected humans 
when different controls are considered. It is possible to see that the upper control, which means 
that everyone is vaccinated, implies that just a few individuals were infected, allowing eradication of 
the disease. Although the optimal control, in the sense of objective ((2]), allows the occurrence of an 
outbreak, the number of infected individuals is much lower when compared with a situation where 
no one is vaccinated. Also, the costs are very low when compared with the upper control case. 
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4 Conclusions 



Dengue is an infectious tropical disease difficult to prevent and manage. Researchers agree that the 
development of a vaccine for dengue is a question of high priority. In the present study we show how 
a vaccine results in saving lives and at the same time in a reduction of the budget related with the 
disease. As future work we intend to study the interaction of a dengue vaccine with other kinds of 
control already investigated in the literature, such as insecticide and educational campaigns [5J |5] . 
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